
Infection Ecology & Epidemiology 


ISSN: (Print) 2000-8686 (Online) Journal homepage: https://www.tandfonline.com/loi/ziee20 



Taylor &. Francis 

Taylor &. Francis Group 


Prevalence and genetic diversity of coronaviruses 
in wild birds, Finland 

Satu Hepojoki, Erika Lindh, Olli Vapalahti & Anita Huovilainen 


To cite this article: Satu Hepojoki, Erika Lindh, Olli Vapalahti & Anita Huovilainen (2017) 
Prevalence and genetic diversity of coronaviruses in wild birds, Finland, Infection Ecology & 
Epidemiology, 7:1, 1408360, DOI: 10.1080/20008686.2017.1408360 

To link to this article: https://doi.Org/10.1080/20008686.2017.1408360 



© 2017 The Author(s). Published by Informa 
UK Limited, trading as Taylor & Francis 
Group. 



View supplementary material GF 



Published online: 28 Nov 2017. 



Submit your article to this journal GF 



Article views: 280 


CrossMark 


View Crossmark data GF 


Full Terms & Conditions of access and use can be found at 
https://www.tandfonline.com/action/joumallnformation?joumalCode=ziee20 



















INFECTION ECOLOGY & EPIDEMIOLOGY, 2017 
VOL 7, 1408360 

https://doi.Org/10.1080/20008686.2017.1408360 



Taylor & Francis 

Taylor &. Francis Group 


RESEARCH ARTICLE 


9 OPEN ACCESS 


0 Check for updates 


Prevalence and genetic diversity of coronaviruses in wild birds, Finland 

Satu Hepojoki 3 * *, Erika Lindh ab *, Olli Vapalahti abc and Anita Huovilainen d 

a Virology, University of Helsinki, Helsinki, Finland; b Helsinki University Hospital, Helsinki, Finland; Veterinary Biosciences, University of 
Helsinki, Helsinki, Finland; d Research and Laboratory Department, Veterinary Virology, Finnish Food Safety Authority, Evira, Helsinki, 
Finland 


ABSTRACT 

Introduction: Migratory birds act as hosts for a number of zoonotic viruses, and have the 
ability to disperse these viruses to distant geographic locations. Coronaviruses (CoVs) repre¬ 
sent a family of zoonotic viruses with wide variety of animal hosts, including birds and 
humans. The infections caused by coronaviruses vary from mild to severe, depending on 
the viral species and the host. Since the coronaviruses exhibit extraordinary large RNA 
genome, also the rate of homologous recombination is high, which in turn contributes to 
the genetic diversity and interspecies host-switches of CoVs. The emergence of novel CoVs 
has been rich during the last decades, and wild birds seem to serve as reservoirs for a variety 
of CoV strains. We examined the CoVs circulating among wild birds in Finland. 

Materials and methods: Samples (cloacal swab, tracheal swab, oropharyngeal swab, or 
tissue) representing 61 bird species were collected during 2010-2013, and examined by RT- 
PCR targeting the RdRp gene for the presence of CoV RNA. 

Results: Altogether 51/939 (5.4%) of the examined birds were found positive by RT-PCR. 
Diverse gamma- and deltacoronavirus sequences were detected. 

Discussion: Gamma- and deltacoronaviruses circulate among wild birds in Finland. The 
number of CoV-positive birds detected each year varies greatly. 
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Introduction 

Coronaviruses (CoVs) comprise a family under the 
order Nidovirales (family Coronaviridae) and infect 
a wide variety of mammals and birds. The course of 
infection varies greatly from asymptomatic to severe 
disease, depending on the host and virus species in 
question. The genome of CoVs is one of the largest 
(25-32 kb) viral RNA-genomes [1]. Based on phy¬ 
logenetic analysis, the CoVs are divided into four 
different genera: Alpha-, Beta-, Gamma-, and 
Deltacoronavirus. The alpha- and betacoronaviruses 
are carried by mammals, whereas the gamma- and 
deltacoronaviruses mainly infect birds, with few 
exceptions [2,3]. The large genomes, infidelity of 
the RNA-dependent RNA polymerase, and high fre¬ 
quency of homologous RNA recombination are the 
main factors contributing to the high genetic diver¬ 
sity of CoVs [4-6]. 

The first CoV, Infectious bronchitis virus (IBV), 
was identified in 1937 [7]. IBV mainly infects chick¬ 
ens, but may infect other bird species as well. IBV is 
highly contagious and affects the respiratory tract, 
gut, kidney, and reproductive systems, causing sub¬ 
stantial economic losses in the poultry industry [8]. 
Despite the global distribution of IBV, poultry in 


Finland remained free of clinical cases until April 
2011 [9] after which outbreaks involving several 
CoV genotypes have occurred in Southern Finland. 

The first human CoVs were identified in 1960s [10- 
12]. The human CoVs cause generally mild to moder¬ 
ate upper respiratory tract infections [13-15]. In 2003, 
a novel highly pathogenic betacoronavirus emerged in 
China, causing severe disease characterized by acute 
respiratory distress and it became known as severe 
acute respiratory syndrome (SARS)-CoV. The emer¬ 
gence of SARS-CoV [16] inspired virologists to more 
explore the highly divergent group of coronaviruses 
and their hosts, leading to the identification of a 
rapidly growing number of CoV species particularly 
in bats [2]. More recently, another highly pathogenic 
betaCoV infecting humans, the Middle East 
Respiratory Syndrome (MERS) CoV emerged in 2012 
with a case fatality rate of over 40% [17,18]. 

Migratory birds have the ability to facilitate the 
dispersion of microorganisms with zoonotic poten¬ 
tial. Wild birds have been associated with the ecology 
and dispersal of at least West Nile virus, tick-borne 
encephalitis virus, influenza A virus (IAV) and 
Newcastle disease virus (NDV) [19-21]. Since the 
discovery of IBV in 1937, it remained the only 
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known Gammacoronavirus for over 50 years, but the 
number has increased dramatically during the last 
10 years [5]. Thereafter, representatives of the genera 
Gamma- and Deltacoronavirus have been isolated 
from both wild and domestic birds including species 
from the order Anseriformes, Pelecaniformes, 
Ciconiiformes, Galliformes, Columbiformes , and 
Charadriiformes [22-24]. In this report we provide a 
description of CoV species circulating in wild birds in 
Finland. Altogether 939 samples representing 61 dif¬ 
ferent bird species were collected during 2010-2013 
and examined for the presence of CoV RNA. 

Materials and methods 
Ethical statement 

All active and passive surveillance samples tested in 
this study, including samples provided by hunters, and 
injured or diseased birds were initially sent to the 
Finnish Food Safety Authority Evira (Helsinki, 
Finland) for avian influenza testing as part of a 
national surveillance program. All hunters had appro¬ 
priate permits, and no birds were killed for research 
purposes. The hunted birds were shot during the 
annual duck-hunting season, from non-endangered 
species, and thereby do not require ethical approval. 


Sampling 

The samples were collected as part of a national 
active and passive surveillance program, coordi¬ 
nated by the Finnish Food Safety Authority Evira 
and University of Helsinki. The sample panels 
included (1) hunted clinically healthy birds (active 
surveillance), (2) hunted clinically diseased birds 
(active surveillance), and (3) birds found dead (pas¬ 
sive surveillance) or (4) clinically diseased (passive 
surveillance) (Table SI). The sample material used 
for RNA detection was cloacal swab, tracheal swab, 
oropharyngeal swab, or tissue. The swab samples 
were collected using nylon swabs and stored and 
transported in Universal Transport Medium (both 
by Copan International). Prior to RNA isolation 
from the swab specimen (stored at -80°C), the 
samples were centrifuged to remove solid particles 
and supplemented with additional antibiotics 
(streptomycin-penicillin). The tissue samples were 
homogenized and separated by centrifugation and 
supplied with antibiotics (streptomycin-penicillin). 
The sample panel from 2010 consisted of 343 sam¬ 
ples, the 2011 panel of 171 samples, the 2012 panel 
of 287 samples, and the 2013 panel of 138 samples. 
The samples were also screened for influenza A 
virus RNA [25]. 


Detection of coronavirus RNA 

QIAamp Viral RNA Mini Kit was used to extract the 
RNA from the samples. The CoV screening was per¬ 
formed using degenerate primers [26] to amplify a 
179 nt region in the polymerase (RdRp) gene (Orflb). 
All positive samples were further characterized by 
amplifying a longer (608-610 bp) region in the 
RdRp gene using previously described primers (for¬ 
ward primer 5'-TGGGWTGGGAYTAYCC 
WAARTGYGA-3' and reverse primer 5'- 
GC AT W GTRT G YT GN G ARC ARA ATT C - 3') [23]. 
All PCR reactions were performed on a BioRad 
PTC-100 Thermal cycler or Applied Biosystems 
Veriti Thermal cycler by a one-step RT-PCR protocol 
using QIAGEN One-Step RT-PCR Kit (Qiagen, 
Germany). The reaction volume for PCR was 25 pi 
containing 5 pi 5x buffer, 1 pi dNTP mix, 1 pM both 
primers, 0.2 pi RNAse inhibitor (Applied Biosystems, 
20 U/pl), 2 pi RNA, and 1 pi Qiagen One-Step RT- 
PCR Enzyme Mix with thermal conditions as follows: 
reverse transcription at 50°C for 30 min, initial PCR 
activation at 95°C for 15 min, followed by denatura- 
tion step at for 30 s, annealing at 50°C for 30 s, and 
extension at 72°C for 30 s. The PCR products were 
analyzed on 2% agarose gels and gel purified using 
QIAquick gel extraction kit (Qiagen). 

Sequencing and phylogenetic analysis 

The RT-PCR amplicons were sequenced in both direc¬ 
tions at the Sequencing Unit of Institute for Molecular 
Medicine Finland, (Helsinki), using BigDye v3.1 
chemistry and run by ABI3730xl DNA Analyzer. The 
sequences were assembled in BioEdit and aligned 
using ClustalW [27-29]. The reference sequences 
were obtained through BLAST search [30]. and from 
GenBank [31]. Phylogenetic analyses were based on 
sequence alignments of 464 nucleotides (position 
100-547 of the RdRp gene, by IBV (AY392086) 
count) within Orflab. Sequence data from 46 of our 
samples was included in the analysis and the rest were 
left out because of inadequate sequence quality. The 
phylogenetic tree was constructed from Maximum- 
likelihood phylogenies and bootstrap values calculated 
by 1000 replicates. Gaps were handled by partial dele¬ 
tion. The evolutionary model GTR + G + I was used 
after evaluation of the best-fit model according to 
Bayesian Information Criterion. The alignments, selec¬ 
tion of evolutionary model and phylogenetic analyses 
were performed within Mega5 [29]. 

Results 

Coronavirus RNA in wild birds 

We screened in total 939 samples from 61 bird spe¬ 
cies (Table SI) during a time period of four years 
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Table 1. Bird samples positive for coronavirus RNA. 


Species 

Clinical status (no.) 

% (no. positive/negative) 

Gamma 

Delta 

Co-infections (no.) 

Anas platyrhynchos 

Hunted clinically healthy (27) and hunted clinically 

21.7% (28/129) 

X 


H9N2 (1) + H3N8 (2) 


diseased (1) 





Anas crecca 

Hunted clinically healthy 

16.4% (9/55) 

X 


H3N8 (1) 

Larus argentatus 

Found dead 

9.6% (5/52) 

X 



Chroicocephalus 

Found dead 

4.2% (1/24) 


X 


ridibundus 

Larus fuscus 

Found dead 

8.3% (1/12) 


X 


Columba sp. 

Hunted clinically diseased (1) and found clinically 

3.6% (2/56) 

X 




diseased (1) 





Cygnus cygnus 

Found dead 

3.6% (1/78) 

X 



Clangula hyemalis 

Found dead 

1.3% (4/8) 

X 



Total 


51/412 





(2010-2013) using a traditional, conserved RT-PCR 
targeting a 179 fragment of the RdRp gene (Orflab) 
of all coronavirus lineages. Altogether 51/939 (5.4%) 
of the examined birds were found to be CoV RNA 
positive (Tables 1 and 2), of which 27 were found 
healthy, and 24 dead or diseased. CoV RNA was 
detected in eight species ( Anas platyrhynchos, Anas 
crecca, Clangula hyemalis, Cygnus cygnus, Larus 
argentatus, Chroicocephalus ridibundus, Larus fuscus, 
and Columpa sp.). CoV RNA was most abundant in 
the samples from 2010 (11%) and 2013 (7.2%), 
whereas during 2011 (0%) and 2012 (0.7%) only few 
of the tested samples were found positive. Influenza 
A virus (IAV) was also detected in four of the CoV 
positive samples (Table 1) [25,32]. The sampling 
locations are illustrated in Figure 1. 


Phylogenetic analyses of the detected avian 
coronaviruses 

The samples positive in the initial RT-PCR screen 
were consequently amplified for phylogenetic pur¬ 
poses by the method described by [23]. The number 
of reference strains for the alignments was unfortu¬ 
nately limited as GenBank deposited sequences vary 
in both location and length as a result of the wide 
array of approaches that have been employed for 
sequencing of coronaviruses. Table 2 shows the 
GenBank accession numbers of the viruses illustrated 
in the phylogenetic tree (Figure 2). 

The viruses detected from wild ducks shared up 
to 99% nucleotide identity with other known 
Co Vs of the Gammacoronavirus genus. The phy- 
logenetically most closely related viruses are iso¬ 
lates from ducks in Russia and USA, China [33], 
Hong Kong [24] and South Korea [33,34]. The 
gammacoronaviruses that were detected in gulls 
in this study shared up to 99% sequence identity 
and clustered with previous isolates from gulls in 
Siberia and Alaska [23]. The viruses isolated from 
pigeons ( Columba sp) (Fin6709 and Finll782) and 
an isolate from a long-tailed duck ( Clangula hye¬ 
malis) shared only up to 89% nucleotide identity 


with the currently published strains, with the 
highest nucleotide identity with duck gammacor¬ 
onaviruses. Two isolates, AvCo V/Larus fuscus/ 
Finland/10059/2013 and AvCoV/ Chroicocephalus 
ridibundus/ Finland/10083/2013, share 83% and 
85% nucleotide identity with published avian and 
mammalian deltacoronaviruses, respectively. They 
clustered in the Deltacoronavirus clade in the ML 
trees, most closely to the White-eye CoV isolated 
in Hong Kong (JQ065044) (Figure 2). None of the 
viruses from wild birds in this study showed close 
relatedness to circulating IBV-strains. 


Discussion 

Several viruses, including zoonotic and economic¬ 
ally significant pathogens, are known to circulate 
among wild birds. Co Vs are a large group of 
viruses infecting mammals and birds, including 
infectious bronchitis virus (IBV), which is a highly 
contagious coronavirus infecting chickens. A vari¬ 
ety of coronaviruses are known to circulate among 
wild birds. We set up this study to determine the 
occurrence of coronavirus infections in wild birds 
in Finland. Altogether 939 samples collected 
between 2010 and 2013 from wild birds were 
included in the study. The prevalence was fairly 
high (up to 11%) although with great annual varia¬ 
tion. The sample panels differed to some extent in 
terms of time and geographic distribution 
(Figure 2) and sampled species, which might partly 
explain the annual variation. For example, majority 
of the samples of 2010 were from duck species 
(Anas sp.) sampled for active surveillance purposes. 
In 2011-2013 on the contrary, most samples were 
derived through passive surveillance and included a 
smaller proportion of ducks. However, no firm 
conclusions can be drawn whether the observed 
difference in prevalence is due to some regular 
pattern, or due to the sampling methods (hetero¬ 
geneity of the sample material). 

Some of the CoV positive samples were also posi¬ 
tive for H3N8 and H9N2 IAV RNA. Similar double 
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Figure 1.A map illustrating the sampling locations. The 
sampling locations (sub-regional units, year 2013) are indi¬ 
cated in the Table 2, and are here illustrated by numbers on a 
map of Finland. 


infections have been reported by other researchers 
[35], and low pathogenic avian influenza viruses 
(LPAI) together with type-1 avian paramyxoviruses 
(APMV) are also commonly detected [36]. The co- 
infections observed here were in four apparently 
healthy ducks, three mallards and a teal, suggesting 
that infections by a number of low virulence viruses 
are well tolerated by migratory ducks. 

CoV infections in mallard ducks ( Anas platyr- 
hynchos) have been documented in Sweden 
[26,35], and one of the studies reported the pre¬ 
valence of CoV infections to be 6.9% albeit with 
seasonal variation [35]. In a more recent study, 
Wille et al. documented the prevalence of Co Vs 
among wild waterbirds in Sweden to be even 
18.7% [37]. A study from Norway reported the 
prevalence of CoV RNA among graylag goose 
(.Anser anser) in Northern Europe to be as high 
as 38% in 2004, but they also documented great 
annual and geographical variation; in 2003 the 
prevalence was only 18% [38]. Also in England, 
wildfowl (Anser if or mes) and waders 
(Charadriiformes) have been reported to carry 
CoVs [39]. A report from South Korea documen¬ 
ted a prevalence of CoV RNA to be 0.95% in 
wildfowl (. Anseriformes ), and they also observed 
annual variation in the detection of CoV RNA 
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Figure 2. Maximum likelihood phylogeny of avian coronaviruses detected in Finland in wild birds. A 464 nucleotide long 
sequence region of the RdRp gene within Orflab was sequenced and aligned with reference sequences: the analysis includes 
alignments of 46 viruses detected in this study and a diverse set of reference strains obtained through BLAST search. Maximum- 
likelihood analysis using 1000 bootstrap replicates was used to infer tree topology. The evolutionary model GTR+G + I was used 
after evaluation of the best-fit model according to Bayesian Information Criterion. Bootstrap support values exceeding 70 are 
shown next to the nodes. Strains described in this study are marked by dots. Branches have been collapsed where our viruses 
share high sequence identity, and the branch wideness is relative to the number of viruses included. 


[34]. Based on the studies from South Korea [34] 
and Northern Europe [38], and the results we 
obtained, it seems that large annual variation in 
the circulating RNA seems to be characteristic to 
avian coronaviruses. A study in China reported 
Gammacoronavirus and Deltacoronavirus RNA in 
12% of the screened, apparently healthy, wild 
aquatic birds ( Anseriformes , Pelecaniformes , and 
Ciconiiformes) [24], which is in concordance with 
our results. Another study, conducted in the 
Beringia area (encompassing areas of Alaska and 
Siberia), reported Gammacoronavirus RNA in 
6.4% of the examined birds (. Anseriformes , 


Pelecaniformes , and Charadriiformes) [23]. 
Interestingly, some of the sequences found in the 
Chinese study [24] were identical to those found 
from the Beringia area [23]. Also the Finnish Duck 
Co Vs reported in this study cluster very closely 
with duck Co Vs of the Gammacoronavirus genus 
from Siberia and China, which are connected by 
migratory routes. The close genetic relatedness of 
these strains (as well as those found in South 
Korea and England) indicates that migratory 
birds do have a specific role in dispersing CoV 
to distinct geographical locations. However, the 
sequence analyses might be somewhat biased due 
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to the limited availability of reference strains. This 
and previous studies have mainly been focused on 
duck and shorebird species of the Northern hemi¬ 
sphere, hence the involvement of other avian spe¬ 
cies and a wider global distribution of avian Co Vs 
remains unclear. 

Interestingly enough our results show that the genus 
Larus (order Charadriiformes) can serve as a host for 
both genera, the Deltacoronavirus (species: Larus fus- 
cus) and Gammacoronavirus (species: Larus argenatus). 
Also, the fact that we detected a Deltacoronavirus in 
Chroicocephalus ridibundus is interesting, since in 
another study [23] Gammacoronavirus RNA was 
found from the same bird species. It was previously 
reported, that the bird species from the order 
Anseriformes can be infected with both delta- and gam- 
macoronaviruses [24], and we now demonstrate that 
the same applies for the order Charadriiformes. The 
abundance of the detected CoV sequences originated 
from ducks supports the hypothesis that the Duck Co Vs 
would persist as endemic in duck populations causing 
subclinical infections similarly to LPAI and APMV. 

To conclude, we show that there is a wide range of 
gammacoronaviruses and also some deltacoronaviruses 
circulating in wild birds in Finland. In addition, we 
show that the number of CoV-positive birds detected 
each year varies greatly. Due to the high recombination 
rate of Co Vs, new species constantly arise that are able 
to adapt to new hosts and ecological niches [5] . Since 
the number of CoV sequences detected from birds was 
so abundant, it would be interesting to study what types 
of Co Vs are circulating in small Finnish mammals, 
especially those in close contact with birds. 

Acknowledgments 

We would like to thank the technical staff members, espe¬ 
cially Merja Hautala, at the Finnish Food Safety Authority 
(Evira) for skillful technical assistance in the laboratory. 
The study was funded by Evira, The Finnish Ministry of 
Agriculture and Forestry MAKERA foundation (grant 
1771/312/2014), and Helsinki University Hospital Funds 
Project (TYH2014251). 

Notes on contributors 

Satu Hepojoki is a postdoctoral researcher based in 
Helsinki. Her main research interests include zoonotic 
viruses, infectious diseases, and rapid diagnostics. 

Erika Lindh has worked in the viral zoonoses research 
group since 2006. She defended her thesis on avian influ¬ 
enza viruses in 2015, at the University of Helsinki, Faculty 
of Medicine. 

Olli Vapalahti is by education MD and PhD and Specialist 
in Clinical Microbiology. He is Professor of Zoonotic 
Virology shared by medical and veterinary faculties of 
University of Helsinki and principal investigator of a 
research group targeting viral zoonoses, and Chief 


Physician at Virology and Immunology, HUSFAB, 
Helsinki University Hospital. 

Anita Huovilainen is a PhD, senior scientist at the 
Virology Department, Finnish Food Safety Authority, 
Helsinki. She also holds the title of virology docent at the 
University of Helsinki. Her main research interest is veter¬ 
inary virology focused on viruses affecting both poultry 
and wild birds. 

Disclosure statement 

No potential conflict of interest was reported by the 
authors. 

Funding 

This work was supported by the Helsinki University 
Hospital Funds Project [TYH2014251]; the Finnish Food 
Safety Authority (Evira); The Finnish Ministry of 
Agriculture and Forestry MAKERA Foundation [1771/ 
312/2014]. 

Notes on contributors 

Satu Hepojoki is a postdoctoral researcher based in 
Helsinki. Her main research interests include zoonotic 
viruses, infectious diseases, and rapid diagnostics. 

Erika Lindh has worked in the viral zoonoses research 
group since 2006. She defended her thesis on avian influ¬ 
enza viruses in 2015, at the University of Helsinki, Faculty 
of Medicine. 

Olli Vapalahti is by education MD and PhD and Specialist 
in Clinical Microbiology. He is Professor of Zoonotic 
Virology shared by medical and veterinary faculties of 
University of Helsinki and principal investigator of a 
research group targeting viral zoonoses, and Chief 
Physician at Virology and Immunology, HUSFAB, 
Helsinki University Hospital. 

Anita Huovilainen is a PhD, senior scientist at the 
Virology Department, Finnish Food Safety Authority, 
Helsinki. She also holds the title of virology docent at the 
University of Helsinki. Her main research interest is veter¬ 
inary virology focused on viruses affecting both poultry 
and wild birds. 


References 

[1] Flint S, Enquist F, Racaniello V, et al., editors. 
Principles of virology. 3rd ed. Washington (USA): 
ASM Press; 2008. 

[2] King AMQ, Adams MJ, Carstens EB, et al., editors. 
Virus taxonomy: classification and nomenclature of 
viruses: ninth report of the international committee 
on taxonomy of viruses. 9th ed. San Diego (CA): 
Academic Press; 2012. 

[3] Berson SA, Yalow RS, Bauman A, et al. Insulin- 
1131 metabolism in human subjects: demonstration 
of insulin binding globulin in the circulation of 
insulin treated subjects. J Clin Invest. 
1956;35:170-190. 


8 @ S. HEPOJOKI ET AL. 


[4] Jackwood MW, Hall D, Handel A. Molecular evolu¬ 
tion and emergence of avian gammacoronaviruses. 
Infect Genet Evol. 2012;12:1305-1311. 

[5] Woo PC, Lau SK, Huang Y, et al. Coronavirus diver¬ 
sity, phylogeny and interspecies jumping. Exp Biol 
Med (Maywood). 2009;234:1117-1127. 

[6] Denison MR, Graham RL, Donaldson EF, et al. 
Coronaviruses: an RNA proofreading machine regu¬ 
lates replication fidelity and diversity. RNA Biol. 
2011;8:270-279. 

[7] Beaudette FR, Hudson CB. Cultivation of the virus of 
infectious bronchitis. J Am Vet Med Assoc. 
1937;90:51-60. 

[8] Cavanagh D. Coronavirus avian infectious bronchitis 
virus. Vet Res. 2007;38:281-297. 

[9] Pohjola LK, Ek-Kommonen SC, Tammiranta NE, 
et al. Emergence of avian infectious bronchitis in a 
non-vaccinating country. Avian Pathol. 2014;43:244- 
248. 

[10] Tyrrell DA, Bynoe ML. Cultivation of viruses from a 
high proportion of patients with colds. Lancet. 
1966;1:76-77. 

[11] Hamre D, Procknow JJ. A new virus isolated from the 
human respiratory tract. Proc Soc Exp Biol Med. 
1966;121:190-193. 

[12] McIntosh K, Dees JH, Becker WB, et al. Recovery in 
tracheal organ cultures of novel viruses from patients 
with respiratory disease. Proc Natl Acad Sci USA. 
1967;57:933-940. 

[13] van der Hoek L, Pyre K, Jebbink MF, et al. 
Identification of a new human coronavirus. Nat 
Med. 2004;10:368-373. 

[14] Woo PC, Lau SK, Chu CM, et al. Characterization and 
complete genome sequence of a novel coronavirus, 
coronavirus HKU1, from patients with pneumonia. J 
Virol. 2005;79:884-895. 

[15] van der Hoek L. Human coronaviruses: what do they 
cause? Antivir Ther. 2007;12:651-658. 

[16] Drosten C, Gunther S, Preiser W, et al. Identification 
of a novel coronavirus in patients with severe acute 
respiratory syndrome. N Engl J Med. 2003;348:1967- 
1976. 

[17] Corman VM, Eckerle I, Bleicker T, et al. Detection of 
a novel human coronavirus by real-time reverse-tran¬ 
scription polymerase chain reaction. Euro Surveill. 
2012;17:20285. 

[18] Cotten M, Watson SJ, Zumla Al, et al. Spread, circula¬ 
tion, and evolution of the middle east respiratory 
syndrome coronavirus. MBio. 2014;5:e01062-13. 

[19] Reed KD, Meece JK, Henkel JS, et al. Birds, migration 
and emerging zoonoses: west nile virus, lyme disease, 
influenza A and enteropathogens. Clin Med Res. 
2003;1:5-12. 

[20] Hubalek Z. An annotated checklist of pathogenic 
microorganisms associated with migratory birds. J 
Wildl Dis. 2004;40:639-659. 

[21] Gilbert M, Xiao X, Domenech J, et al. Anatidae migra¬ 
tion in the western Palearctic and spread of highly 
pathogenic avian influenza H5NI virus. Emerg Infect 
Dis. 2006;12:1650-1656. 

[22] Honkavuori KS, Briese T, Krauss S, et al. Novel cor¬ 
onavirus and astrovirus in Delaware Bay shorebirds. 
PLoS One. 2014;9:e93395. 


[23] Muradrasoli S, Balint A, Wahlgren J, et al. Prevalence 
and phylogeny of coronaviruses in wild birds from the 
Bering Strait area (Beringia). PLoS One. 2010;5: 
el3640. 

[24] Chu DK, Leung CY, Gilbert M, et al. Avian corona¬ 
virus in wild aquatic birds. J Virol. 201 1;85:12815— 
12820. 

[25] Lindh E. Avian influenza and newcastle disease 
viruses in Finland - genetics, epidemiology and ecol¬ 
ogy in the natural host, wild waterfowl [Ph.D. thesis]. 
Helsinki: Faculty of Medicine, University of Helsinki; 
2015. 

[26] Muradrasoli S, Mohamed N, Hornyak A, et al. 
Broadly targeted multiprobe QPCR for detection of 
coronaviruses: coronavirus is common among mallard 
ducks (Anas platyrhynchos). J Virol Methods. 
2009;159:277-287. 

[27] Larkin MA, Blackshields G, Brown NP, et al. Clustal 
W and Clustal X version 2.0. Bioinformatics. 
2007;23:2947-2948. 

[28] Hall TA. BioEdit: a user-friendly biological 
sequence alignment editor and analysis program 
for Windows 95/98/NT. Nucl Acids Symp Ser. 
1999;41:95-98. 

[29] Tamura K, Stecher G, Peterson D, et al. MEGA6: 
molecular evolutionary genetics analysis version 6.0. 
Mol Biol Evol. 2013;30:2725-2729. 

[30] The Basic Local Alignment Search Tool (BLAST), 
[cited 2016 Jul 3]. Available from: http://blast.ncbi. 
nlm. nih. go v/Blast. cgi 

[31] GenBank NCBI. [cited 2016 Jul 3] Available from: 
http://www.ncbi.nlm.nih.gov/genbank/ 

[32] Lindh E, Ek-Kommonen C, Vaananen V-M, et al. 
Molecular epidemiology of H9N2 influenza viruses in 
Northern Europe. Vet Microbiol. 2014;172:548-554. 

[33] Zhuang Q-Y, Wang K-C, Liu S, et al. Genomic ana¬ 
lysis and surveillance of the coronavirus dominant in 
ducks in China. PLoS One. 2015;10:e0129256. 

[34] Kim H-R, Oem J-K. Surveillance of avian corona¬ 
viruses in wild bird populations of Korea. J Wildl 
Dis. 2014;50:964-968. 

[35] Wille M, Avril A, Tolf C, et al. Temporal dynamics, 
diversity, and interplay in three components of the 
virodiversity of a Mallard population: influenza A 
virus, avian paramyxovirus and avian coronavirus. 
Infect Genet Evol. 2015;29:129-137. 

[36] Tolf C, Wille M, Haidar A-K, et al. Prevalence of avian 
paramyxovirus type 1 in Mallards during autumn 
migration in the western Baltic Sea region. Virol J. 
2013; 10:285-422X-10-285. 

[37] Wille M, Muradrasoli S, Nilsson A, et al. High pre¬ 
valence and putative lineage maintenance of avian 
coronaviruses in Scandinavian waterfowl. PLoS One. 
2016;ll:e0150198. 

[38] Jonassen CM, Kofstad T, Larsen I-L, et al. Molecular 
identification and characterization of novel corona¬ 
viruses infecting graylag geese (Anser anser), feral 
pigeons (Columbia livia) and mallards (Anas platyr¬ 
hynchos). J Gen Virol. 2005;86:1597-1607. 

[39] Hughes LA, Savage C, Naylor C, et al. Genetically 
diverse coronaviruses in wild bird populations of 
northern England. Emerg Infect Dis. 2009; 15:1091- 
1094. 


